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The  lattice  spacing  comparator  estab- 
lished at  the  National  Institute  of  Stan- 
dards and  Technology  to  measure  the 
lattice  spacing  differences  between 
nearly  perfect  crystals  is  described  in 
detail.  Lattice  spacing  differences  are 
inferred  from  the  measured  differences 
in  Bragg  angles  for  different  crystals. 
The  comparator  is  a two  crystal  spec- 
trometer used  in  the  nearly  nondisper- 
sive  geometry.  It  has  two  x-ray  sources, 
two  detectors,  and  a device  which  per- 
mits remote  interchange  of  the  second 
crystal  sample.  A sensitive  heterodyne 
interferometer  which  is  calibrated  with 
an  optical  polygon  is  used  to  measure 
the  Bragg  angles.  The  crystals  are  man- 
ufactured with  nearly  equal  thicknesses 
so  that  the  recorded  profiles  exhibit 
pendellosung  oscillations  which  permit 


more  precise  division  of  the  x-ray 
profiles.  The  difference  in  lattice  spac- 
ing between  silicon  samples  used  at 
Physikalisch-Technische  Bundesanstalt 
(PTB)  and  NIST  has  been  measured 
with  a relative  uncertainty  of  IxlO'**. 
This  measurement  is  consistent  with  ab- 
solute lattice  spacing  measurements 
made  at  PTB  and  NIST.  Components 
of  uncertainty  associated  with  system- 
atic effects  due  to  misalignments  are 
derived  and  estimated. 
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1.  Introduction 

Absolute  measurements  of  the  crystal  lattice 
spacing  of  silicon  have  been  reported  by  two  labo- 
ratories [1, 2,3,4]  and  are  in  progress  in  at  least  two 
other  laboratories  [5,6].  These  experiments  present 
a very  difficult  metrological  challenge  and  require 
extreme  care  in  order  to  achieve  a relative  uncer- 
tainty of  several  parts  in  10**.  The  output  of  these 
experiments  is  the  lattice  spacing  of  a particular 
silicon  specimen  in  terms  of  an  optical  wavelength 
standard. 
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Crystals  having  well  defined  lattice  spacings  and 
whose  lattice  spacings  are  known  with  an  uncer- 
tainty in  the  0.1  to  0.01  ppm  range  (lppm  = 10"'’) 
are  needed  in  the  field  of  metrology  (i.e.,  x-ray  and 
gamma-ray  wavelength  measurements).  However, 
the'effort  required  to  directly  compare  a large  fam- 
ily of  crystal  samples  to  an  optical  wavelength  stan- 
dard is  prohibitive.  Thus,  several  techniques  and 
devices  to  more  rapidly  measure  lattice  spacing  dif- 
ferences with  a relative  uncertainty  below  0.1  ppm 
have  been  devised.  This  paper  describes  in  detail 
the  precision  lattice  spacing  comparator  (At/  com- 
parator) established  at  the  National  Institute  of 
Standards  and  Technology  and  reports  the  results 
of  a careful  comparison  of  Si  samples  which  were 
measured  absolutely  by  the  x-ray  interferometer 


technique  at  the  Physikalisch-Technische  Bunde- 
sanstalt  (PTB)  and  NIST.  The  measured  lattice 
spacing  difference  has  a relative  uncertainty  of 
1.4  X 10"*  and  is,  at  present,  consistent  with  the  ab- 
solute measurements  reported  by  the  two  laborato- 
ries. 

2.  Principle  of  the  Measurement 

The  principle  of  the  NIST  lattice-comparison 
measurement  is  illustrated  by  the  geometrical  lay- 
out of  the  spectrometer  (Fig.  1).  This  geometry, 
which  is  similar  to  one  suggested  by  Hart  [7],  is  a 
conventional  double-crystal  Laue-case  spectrome- 
ter. Ando  et  al.  [8],  Becker  et  al.  [9],  and  Hauser- 
mann  et  al.  [10]  have  also  used  similar  geometries 


Fig.  1.  Schematic  diagram  of  the  NIST  lattice  spacing  comparator.  Si,  S2  — x-ray  sources;  Li,  L2  — shutters;  Aj,  A2  — apertures;  Ci, 
C2  — crystals;  Di,  D2  — detectors;  /,  /’  Ar-axis  projections  of  the  x-ray  paths  from  Aj  to  Ci,  Ci  to  C2,  C2  to  A2,  respectively,  a)  Top  view 
{x,  y plane),  b)  side  view  {x,  z plane) 
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to  precisely  compare  lattice  spacings.  As  will 
emerge  in  the  discussion  which  follows,  the  NIST 
instrument  has  certain  new  features  which  enhance 
its  generality  and  sensitivity.  The  spectrometer  has 
two  sources,  Si  and  S2,  two  crystals,  Ci  and  C2  and 
two  detectors,  Di  and  D2,  and  is  used  to  record 
only  nearly  nondispersive  profiles.  The  x-ray  beams 
are  defined  by  precision  source,  Ai,  and  detector, 
A2,  apertures  and  there  are  computer  controlled 
shutters,  Li  and  L2,  near  the  sources  and  the  detec- 
tors, respectively. 

By  using  the  spectrometer  in  the  nearly  nondis- 
persive mode,  the  large  intrinsic  linewidths  associ- 
ated with  characteristic  x-ray  emission  do  not 
contribute  appreciably  to  the  width  of  the  recorded 
profiles.  The  source  shutters  permit  isolation  of  the 
two  x-ray  paths  so  that  doubly  diffracted  profiles 
from  one  source  can  be  recorded  without  interfer- 
ence from  singly  diffracted  profiles  from  the  other 
source.  There  are  three  apertures  in  front  of  each 
detector  (Fig.  lb),  one  above,  one  on,  and  one  be- 
low the  plane  of  dispersion.  An  indexed  shutter  se- 
lects one  of  these  apertures  depending  on  whether 
lattice  spacing  difference  scans  or  crystal  alignment 
scans  are  being  recorded.  For  lattice  spacing  differ- 
ence scans,  apertures  on  the  plane  of  dispersion 
are  used  for  both  sources  and  the  detectors.  As 
crystal  1 is  rotated  through  the  reflection,  the 
source  shutters  alternately  pass  and  block  the  two 

x-ray  paths,  path  (•  • •)  and  path  ( ),  so  that  the 

two  x-ray  profiles  are  essentially  recorded  simulta- 
neously. Thus  errors  due  to  drifts  of  the  crystals  or 
the  angle  interferometer  are  eliminated.  For  crystal 
alignment  scans,  only  one  source  is  used  and  the 
detector  shutter  alternately  permits  radiation  to 
pass  through  the  apertures  above  and  below  the 
plane  of  dispersion  as  crystal  1 is  scanned.  The  an- 
gular shift  of  these  two  out-of-plane  profiles  is  re- 
lated to  the  tilt  misalignment  of  the  planes  of  the 
first  crystal  with  respect  to  the  planes  of  the  second 
crystal.  By  careful  alignment,  the  contribution  to 
the  uncertainty  due  to  crystal  tilt  misalignment  can 
be  made  small  (see  Appendix  A). 

General  features  of  the  comparator  operation 
may  be  understood  as  follows:  If  the  crystal  planes 
near  the  ends  of  the  first  crystal  are  parallel  and 
the  lattice  spacings  of  the  first  (di)  and  the  second 
{d2)  crystals  are  equal  (^1  = ^2),  then  the  profiles  in 

the  path  (•  • •)  and  the  path  ( ) will  peak  at  the 

same  angular  setting  of  the  first  crystal.  However,  if 
di^d2,  then  there  will  be  an  angular  offset,  /32,i  be- 
tween the  two  profiles  which  is  a measure  of  the 
lattice  spacing  difference,  =d2-di,  between  the 
two  crystals.  If  9i  and  62  are  the  Bragg  angles  of 


crystals  1 and  2,  respectively,  then  ^.\  = 2{&2  - 6i)  = 
2A0.  If  Ad/d  is  small  (<10“^)  and  the  required 
uncertainty  in  Ad/d  is  ~ 10“^  then  it  is  appropriate 
to  use  the  differential  form  of  the  Bragg  equation 
to  express  Ad/d. 

Ad/d  = -(AO)  cote  = -{(32.i/2)  cote  (1) 

In  our  comparison  scheme,  the  first  long  crystal 
serves  as  a temporary  (a  few  hours)  reference.  Pro- 
files are  recorded  using  one  of  the  small  second 
crystals.  Then  another  small  second  crystal  is 
brought  into  the  position  where  the  x-ray  beams 
diffracted  by  the  first  crystal  intersect.  The  first 
crystal  is  again  scanned  and  profiles  are  recorded. 

By  subtracting  the  angular  offsets,  /3,  measured 
with  the  different  second  crystal  samples,  the  prop- 
erties of  the  first  crystal  such  as  the  absolute  lattice 
spacing,  the  effects  of  variation  in  lattice  spacing, 
and  parallelism  of  crystal  planes  at  the  two  ends  of 
the  crystal  are  eliminated.  One  needs  only  assume 
that  the  properties  of  the  first  crystal  are  constant 
during  the  time  of  measurement  of  the  two  differ- 
ent second  crystal  samples.  The  equation  for  the 
difference  in  crystal  lattice  spacings  for  the  pair  of 
small  crystals  becomes 

Ad/ds=  - (AO)  cot  Os  = -(13/2)  cotS^ , (2) 

where  ft  = the  Bragg  angle  for  the  standard  sample, 
ft==the  Bragg  angle  for  the  unknown  sample,  ds  = 
lattice  spacing  of  standard  sample,  d^  = lattice  spac- 
ing of  unknown  sample.  Ad  =ds-du,  /3s  = 2(ft-ft) 
= angular  offset  for  the  standard  sample, 
/3u  = 2(0u- ft)  = angular  offset  for  the  unknown 
sample,  and  )3  = )3s  — /3u  = 2( ft  - ft)  = 2A  0. 

In  Appends  A the  equations  for  calculating  Ad 
in  ideal  and  misaligned  geometries  are  derived.  For 
small  lattice  parameter  differences  and  well 
aligned  x-ray  beams  and  crystals,  the  more  involved 
equations  derived  in  the  appendix  reduce  to  Eq. 
(2).  The  forms  of  the  corrections  associated  with 
misalignments  are  explicitly  given  and  their  magni- 
tudes are  estimated  for  our  particular  spectrome- 
ter. This  appendix  is  likely  to  be  most  interesting  to 
the  precision  x-ray  specialist.  Throughout  this  pa- 
per and  particularly  in  the  appendix  the  path  (•  • •) 

is  referred  to  as  the  path  (-)  and  the  path  ( ) 

is  referred  to  as  the  path  ( -I- ). 

Since  x rays  from  the  two  paths  interrogate  the 
same  area  of  the  second  crystal,  small  variations  in 
lattice  spacing  along  a crystal  can  be  measured  with 
this  scheme.  The  second  crystal  samples  reside  on 
a precision  slide  which  allows  easy  interchange  of 
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crystal  samples.  A typical  data  sequence  involves 
many  changes  of  crystal  samples  and  permits  com- 
pensation for  drifts  of  the  crystals  and  the  angle 
measuring  spectrometer. 

3.  Experimental  Apparatus 

The  lattice  comparison  spectrometer  rests  on  a 
1 m X 2 m cast  iron  plate  which  is  isolated  from 
building  vibrations  by  three  passive  air  bags.  The 
x-ray  tubes,  detectors,  and  shutters  are  rigidly  at- 
tached to  this  plate.  A smaller  plate  (61  cmx  122 
cm)  on  which  the  laser,  receiver,  and  precision 
spectrometer  are  mounted  rests  on  vibration  damp- 
ing rubber  feet  on  top  of  the  large  cast  iron  plate. 
The  precision  spectrometer  is  constructed  on  a 35 
cm  X 74  cm  X 2.5  cm  thick  cast  iron  U channel. 
The  angle  interferometer,  rotary  table,  and  transla- 
tion table  reside  on  top  of  the  inverted  U channel 
and  the  drive  arm  coupled  to  the  rotary  table  is 
conveniently  located  underneath  the  channel.  The 
important  dimensions  which  define  the  scale  of  the 
spectrometer  are:  1)  x-ray  focal  spot  to  source  slit 
= 6 cm;  2)  source  slit  to  x-ray  tube  pivot  = 25.4  cm; 
3)  x-ray  tube  pivot  to  first  crystal  = 8.25  cm;  4)  first 
crystal  to  second  crystal  = 8.25  cm;  and  5)  second 
crystal  to  detector  slit  = 28.5  cm. 

The  radiation  sources  are  identical  silver  or 
molybdenum  x-ray  tubes  operated  at  K = 40  keV 
and  7 = 15  mA  and  having  a 1 mm  x 1 mm  focal 
spot  (in  projection).  The  tubes  are  cooled  by  tem- 
perature regulated  water  so  that  the  spectrometer 
temperature  remains  constant  and  uniform.  The  x- 
ray  sources  pivot  about  the  point  of  intersection  of 
the  x-ray  beams  which  lies  between  the  sources  and 
the  first  crystal.  Bragg  angles  from  10°  to  22°  are 
accessible.  The  vertically  defining  slits  in  front  of 
the  sources  have  a 1 mm  opening.  The  center  of 
these  slits  and  the  middle  slits  near  the  detectors 
are  carefully  placed  on  the  plane  of  dispersion 
(plane  perpendicular  to  the  axis  of  rotation)  to 
within  0.1  mm.  The  x-ray  tubes  are  positioned  so 
that  the  focal  spots  fill  the  source  and  detector  slits 
uniformly.  The  lead  shutters  near  the  source  slits 
are  computer  controlled  to  block  or  pass  the  x-ray 
beams. 

The  axis  of  rotation  is  a precision  commercial 
rotary  table  [11]  having  a radial  concentricity  and 
axial  movement  less  than  one  micrometer.  The 
table  is  driven  by  a stepping  motor  coupled  to  a 
threadless  screw  and  a 48  cm  long  tangent  arm  with 
fine  adjustment  of  the  rotation  obtained  by 
piezoelectric  transducers.  The  long  first  crystal  is 


mounted  in  the  center  of  the  rotating  table  so  that 
its  diffracting  planes  are  parallel  to  the  axis  of  rota- 
tion to  within  a few  seconds  (of  plane  angle). 

Angles  are  measured  by  a polarization  encoded 
Michelson  interferometer  having  an  angular  sensi- 
tivity of  a few  X 10“'*  seconds.  The  interferometer 
is  shown  schematically  in  Fig.  2 where  the  inset 
shows  the  over  and  under  arrangement  of  the 
beams.  The  interferometer  is  illuminated  by  a com- 
mercial 633  nm  HeNe  laser  [12]  which  emits  two 
orthogonal  linearly  polarized  frequencies  sepa- 
rated by  1.8  MHz.  The  polarization  sensitive  beam 
splitter,  b,  divides  the  incoming  light  beam  so  that 
one  polarization  and  frequency  traverses  one  arm 
of  the  interferometer  and  the  other  polarization 
and  frequency  traverses  the  other  arm.  The  beam 
splitter,  steering  elements  (e,  g,  h),  and  compensat- 
ing element  (f)  are  chosen  and  positioned  so  that 
the  two  interferometer  arms  have  equal  air  and 
equal  glass  paths  at  zero  angle.  Elements  c and  d 
are  90°  polarization  rotators  which  rotate  the  plane 
of  polarization  of  the  outgoing  and  return  beams 
respectively  in  order  that  both  corner  cubes  are 
transversed  with  the  same  polarization  orientation 
and  that  the  output  beams  are  directed  away  from 
the  laser.  The  roof  prisms,  i and  j,  return  the  light 
to  the  corner  cubes,  k and  1,  doubling  the  angular 
sensitivity  of  the  interferometer  and  keeping  the 
output  beams  fbced  in  space.  The  corner  cubes  are 
mounted  on  an  arm  which  is  rigidly  attached  to  the 
rotating  crystal  table  and  retroreflect  the  light  in 
each  arm  of  the  interferometer. 


0.0  0.2  0.4  0.6  0.8  1.0 


crystal  thickness  • mm 

Fig.  2.  Schematic  diagram  of  the  angle  interferometer,  a,  e,  g, 
h — beam  steering  elements;  b — polarization  sensitive  beam 
splitter;  c,  d — 90°  polarization  rotators;  f— glass  path  compen- 
sating plate;  i,  j — roof  prisms;  k,  I — corner  cubes.  The  insert 
shows  the  passage  of  the  beam  through  the  corner  cube  — roof 
prism  part  of  the  interferometer. 
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The  interferometer  output  beam  which  includes 
the  two  orthogonally  polarized  frequencies  is  ana- 
lyzed by  a 45°  polarizer  included  in  a commercial 
detector  which  provides  a measurement  difference 
frequency  signal.  The  laser  and  its  electronics  also 
provides  a reference  difference  frequency  signal 
for  the  light  emitted  by  the  laser.  For  a stationary 
interferometer  the  two  difference  frequency  signals 
are  identical  except  for  a phase  difference  <f),  which 
is  proportional  to  the  difference  in  the  optical  path 
lengths  of  the  two  interferometer  arms.  Rotation  of 
the  interferometer  arm  causes  a phase  shift  of  the 
measurement  difference  frequency  signal  relative 
to  the  reference  difference  frequency  signal.  Con- 
versely, phase  comparison  of  the  measurement  dif- 
ference frequency  signal  with  the  reference 
difference  frequency  signal  permits  measurement 
of  the  angular  rotation. 

The  1.8  MHz  difference  frequency  signals  must 
be  converted  into  signals  which  permit  up-down 
counting  for  clockwise  and  counter  clockwise 
rotation  of  the  axis,  measurement  of  fringe  frac- 
tions, and  servo  control  of  the  axis  position. 
Although  electronics  to  perform  these  functions  is 
commercially  available,  no  existing  system  could  be 
conveniently  coupled  to  our  spectrometer.  An  elec- 
tronics module  to  perform  the  above  functions  has 
been  developed  and  is  briefly  described  in 
Appendbc  B. 

The  diffraction  profiles  are  recorded  as  x-ray  in- 
tensity vs  interferometer  fringes  and,  thus,  the  an- 
gular offsets,  j8  are  also  measured  in  interferometer 
fringes.  However,  the  angular  offsets,  )3,  needed  in 
Eq.  (2)  must  be  in  absolute  angles.  In  order  to  con- 
vert interferometer  fringes  into  absolute  angles,  the 
angle  interferometer  is  calibrated  using  an  optical 
polygon  and  a photoelectric  autocollimator  [4].  A 
24-sided  optical  polygon  (external  angles  —15°)  is 
mounted  on  top  of  an  indexing  table  on  the  rota- 
tion axis  in  place  of  the  first  crystal.  The  photoelec- 
tric autocollimator  senses  the  directions  normal  to 
the  faces  of  the  optical  polygon.  The  24  external 
angles  are  measured  in  terms  of  interferometer 
fringes  and  the  sum  of  the  angles  is  constrained  to 
be  360°.  The  equation  relating  angles  and  fringes  is 
F =Ksin(6)  where  F is  the  fringe  number,  6 is  the 
corresponding  angle  of  rotation,  and  K is  the  cali- 
bration constant.  The  measured  value  of  K is 
5138551.7  at  22.5  °C.  Near  0 = 0°,  the  value  for  the 
angular  rotation/fringe  is  0.0401  arcsec/fringe.  The 
calibration  constant  is  directly  related  to  the  sepa- 
ration of  the  corner  cubes.  Because  the  arm  which 
defines  the  separation  of  the  corner  cubes  is  made 
out  of  stainless  steel,  the  calibration  constant  has  a 


large  temperature  coefficient  ( — 16.7  ppm/°C).  As 
long  as  the  angular  offsets  are  small  and  the  tem- 
perature is  well  controlled,  this  large  temperature 
coefficient  can  be  tolerated.  However,  in  the  near 
future  the  spectrometer  will  be  equipped  with  an 
invar  corner  cube  arm. 

The  zero  angle  of  the  interferometer  is  deter- 
mined by  measuring  one  of  the  external  angles  of 
the  polygon  in  a symmetric  and  asymmetric  fash- 
ion. In  the  symmetric  measurement,  the  polygon  is 
positioned  on  the  axis  so  that  the  faces  forming  the 
15°  angle  are  measured  at  +7.5°  and  -7.5°.  This 
symmetric  measurement  of  the  external  angle  is 
quite  insensitive  to  the  zero  angle  of  the  interfer- 
ometer. In  the  asymmetric  measurement,  the 
polygon  is  positioned  so  that  the  faces  of  the  same 
angle  are  measured  at  ±14.5°  and  +0.5°.  This 
asymmetric  measurement  of  the  external  angle  is 
very  sensitive  to  the  zero  angle  of  the  interferome- 
ter. The  zero  angle  is  determined  by  requiring 
equality  of  the  two  measurements  of  the  same  an- 
gle. The  uncertainty  of  the  zero  angle  deter- 
mination is  approximately  2 arcsec  which  makes  a 
completely  negligible  contribution  to  the  measure- 
ment of  the  small  angular  offsets  near  zero  angle. 

The  second  crystals  reside  on  a precision  transla- 
tor equipped  with  a stepping  motor  and  a linear 
encoder  [13].  The  slide  has  15  cm  of  travel  and  a 
positioning  accuracy  of  0.01  mm.  Although  the 
slide  has  pitch  and  yaw  of  a few  seconds  in  a few 
cm  of  travel,  it  is  reproducible  at  the  0.01  second 
level  and  shows  no  short  term  drift.  Each  second 
crystal  is  mounted  on  a flexure  hinge  so  its  planes 
can  be  oriented  to  be  roughly  parallel  to  the  planes 
of  the  first  crystal.  Fine  adjustment  is  achieved  by 
applying  voltage  to  a piezoelectric  “tipper.” 

The  diffracted  x-rays  are  detected  by  two  identi- 
cal Nal(Tl)  detectors  which  pivot  about  the  second 
crystal  position.  In  front  of  each  detector  are  three 
slits  each  1.6  mm  x 1.6  mm.  The  slits  are  arranged 
in  the  direction  normal  to  the  plane  of  dispersion 
with  center  of  the  middle  opening  on  the  plane  of 
dispersion  and  the  top  and  bottom  openings  3.2 
mm  above  and  below  the  plane  of  dispersion.  Each 
slit  has  a computer  controlled  shutter  which  per- 
mits radiation  to  pass  through  any  one  of  the  three 
openings. 

4.  Crystal  Preparation 

In  order  to  make  measurements  with  an  uncer- 
tainty of  0.01  ppm,  considerable  care  must  be 
exercised  in  cutting,  polishing,  mounting,  and 
aligning  the  crystal  samples.  The  full  width  at  half 
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maximum  (FWHM)  of  nondispersive  x-ray  profiles 
obtained  with  Ag  Ka  radiation  and  the  Si  440  re- 
flection is  about  0.6  seconds.  Since  the  angular  sep- 
aration which  corresponds  to  a lattice  spacing 
difference  of  1 x 10“*  is  ~ 0.001  seconds,  it  appears 
that  peak  positions  must  be  determined  to  1/500  of 
the  FWHM.  This  is  a formidable  task  requiring 
very  reproducible  profiles.  However,  by  carefully 
matching  the  thickness  of  the  first  and  second  crys- 
tals to  within  a few  micrometers,  nondispersive 
Laue-case  rocking  curves  exhibit  oscillatory  fine 
structure  which  has  a modulation  period  which  is 
typically  less  than  1/10  of  the  FWHM  of  the  profile. 
This  fine  structure  is  predicted  by  the  dynamical 
theory  of  x-ray  diffraction  and  has  been  used  in 
high  precision  determinations  of  the  structure  fac- 
tor of  Si  [14,  15]  and  Ge  [16]  and  in  x-ray  refractive 
index  measurents  [17].  We  use  this  fine  structure 
on  the  x-ray  profiles  as  a sharp  convenient  refer- 
ence to  measure  the  angular  separation  of  the  pro- 
files. 

The  integrated  reflectivity  for  the  double  crystal 
profiles  oscillates  as  a function  of  crystal  thickness. 
In  Fig.  3,  the  integrated  reflectivity  vs  crystal  thick- 
ness is  shown  for  the  Si  440  diffraction  of  Ag  Ka 
radiation.  We  have  chosen  a crystal  thickness  of 
0.455  mm  as  a compromise  of  intensity  and  practi- 
cal problems  with  thin  crystals.  Theoretical  and 
experimental  profiles  for  this  thickness  are  shown 
in  Fig.  4. 

The  theoretical  profile  [18,  19]  is  obtained  by 
adding  the  convolution  of  two  intrinsic  reflections 
for  the  two  polarization  states  (cr  and  tt) 

I(A)  = lo'Z  + d0 

a.TT  ^ 


20000 


Fig.  4.  Theoretical  (lower)  and  experimental  (upper)  profiles 
for  the  440  reflection  of  Ag  Ka  radiation  using  0.455  mm  thick 
Si  crystals. 

where  A = the  angle  between  the  diffracting 
planes  of  crystal  1 and  crystal  2 
=exp  ( - ju//cos0B)|sin[/l (y^-l-v^)*^^]/ 

A = Trf/fo 

Xh  — Xhr  -f"  iXhi 

y =zlsin2ftj/(C|;<:hr|) 

= XhXfi/\xhXfi\ 

to  =Acos0B/(C|?i;hr|) 

^ _ I 1 <T  polarization 
1 1 cos20b  I TT  polarization 

and 

0B  is  the  Bragg  angle 
is  the  linear  absorption  coefficient 
t is  the  crystal  thickness 
to  is  the  extinction  length 
y is  the  angular  parameter 
C is  the  polarization  factor 
A is  the  wavelength 

Xh  is  the  complex  electronic  susceptibility  of  the 
crystal:  xor  + ixu 

For  low  absorption  where  xor>Xhi,  it  can  be 
shown  that 

v^  = 1 -i(e^i/o)/(7rcos0B) 
to  =7rFcos0B/(AAeF|C|) 

/’"■'"(  0)  = [exp(  - /i//cos0B)/(>'^+ 1)] 


Fig.  3.  Integrated  reflectivity  of  double  crystal  nondispersive 
profiles  as  a function  of  crystal  thickness  for  the  Si  440  reflection 
and  Ag  Kaj  radiation. 
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where  e = (|xhi|/|?foi|)|C | 

Xhi,  = imaginary  parts  of  the  Fourier  coeffi- 
cients Xh,  Xo 

K=  volume  of  the  unit  cell 

re  = classical  electron  radius 

F = temperature  modified  structure  factor 

When  generating  a theoretical  profile,  numerical 
values  must  be  assigned  to  the  following  quantities 
fi,  t,  F,  A,  Te,  V,  and  e.  Conversely,  when  fitting  an 
experimental  profile  with  the  theoretical  descrip- 
tion, a quantity  such  as  r or  F can  be  varied  until 
the  difference  between  the  experimental  and  theo- 
retical profiles  is  a minimum. 

In  order  to  obtain  strain  free  samples,  the  crys- 
tals are  made  with  a 1.5  cm  x 1.5  cm  base  on  top  of 
which  is  the  thin  1.5  cm  high  wafer  which  is  used 
for  diffraction.  Drawings  of  the  long  first  crystal 
and  the  second  crystal  are  shown  in  Figs.  5 and  6. 
Typical  lengths  of  the  long  and  short  crystals  are 
6.3  cm  and  1.2  cm,  respectively.  The  long  crystal  is 
balanced  on  a 1.5  mm  wide  silicon  rod  which  is 
centered  on  an  optically  polished  silicon  cylinder. 
This  assembly  is  waxed  together  with  a low  temper- 
ature optician’s  wax.  The  polished  silicon  cylinder 
is  attached  to  a ground  cast  iron  base  with  epoxy 
resin.  The  cast  iron  base  and  the  rotating  table  are 
joined  by  screw  fasteners.  After  assembly  the  angu- 
lar offset  between  the  normal  to  the  polished  sili- 
con cylinder  and  the  crystal  planes  is  measured 
using  an  x-ray  spectrometer  and  an  autocollimator. 
The  long  crystal  is  mounted  on  a flexure  hinge  on 
the  rotating  table  and  is  tipped  an  amount  equal  to 
the  measured  offset.  In  this  way  the  crystal  planes 
are  made  parallel  to  the  axis  of  rotation  with  an 
estimated  uncertainty  of  2 seconds. 

The  short  crystal  assembly  consists  of  a ground 
cast  iron  base,  a piezoelectric  tipper,  a silicon  anvil 
(with  a small  hole  for  a thermistor),  and  the  crystal. 
The  base,  PZT  elements,  and  silicon  anvil  are  at- 
tached by  conductive  epoxy  and  the  crystal  is  bal- 
anced on  the  raised  portion  of  the  silicon  anvil  (1.5 
mm  wide)  and  attached  by  low  temperature  opti- 
cian’s wax.  This  assembly  is  mounted  on  a flexure 
hinge  on  the  translation  table.  Using  x rays  dif- 
fracted by  the  first  crystal  above  and  below  the 
plane  of  dispersion,  the  planes  of  the  second  crys- 
tal are  made  parallel  to  those  of  the  first  crystal  by 
the  coarse  mechanical  adjustment  and  the  fine 
PZT  tipper. 

The  crystals  are  cut  with  a diamond  saw  and 
then  polished  to  the  desired  thickness  using  a 
chemically  assisted  mechanical  polishing  solution 
[20].  Typical  polishing  rates  are  10  ixm  per  h and 


50  p,m  to  70  |xm  are  removed  from  each  side  to 
insure  that  no  saw  damage  remains.  Crystal  thick- 
ness was  measured  using  a coordinate  measuring 
machine  [21].  The  reproducibility  of  the  thickness 
measuring  procedure  was  1 |jLm-2  jxm.  Polishing 
the  long  first  crystal  so  that  the  two  areas  used  for 
diffraction  had  equal  thickness  was  the  most  diffi- 
cult task.  Because  the  edges  of  the  crystals  tend  to 
be  too  thin,  the  areas  used  for  diffraction  were  at 
least  2 mm  from  any  crystal  edge.  After  polishing, 
all  surfaces  of  the  crystal  except  for  the  polished 
surfaces  were  etched  in  HF/HNO3  to  relieve  strain. 


Fig.  5.  Detailed  drawings  of  the  first  crystal,  a — cast  iron,  b — 
Si  disk,  and  c — Si  rod,  and  d-Si  diffraction  crystal. 


Fig.  6.  Detailed  drawing  of  the  second  crystal,  a — cast  iron, 
b — PZT,  c — Si  support,  d — Si  diffraction  crystal. 
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5.  Temperature  Measurement  and  Con- 
trol 

Because  the  lattice  parameter  measurements  are 
relative  measurements,  we  required  only  accurate 
relative  temperature  values  and  not  accurate  abso- 
lute values.  Since  the  expansion  coefficient  for  sili- 
con is  2.56x  10"‘'/K,  relative  temperature 
measurements  accurate  to  4 mK  suffice  for  1 x 10“** 
lattice  parameter  measurements.  Six  ultra-stable 
calibrated  thermistor  probes  are  used  to  measure 
the  temperature  of  the  crystals:  two  thermistors  on 
the  base  of  the  long  first  crystal  (one  on  each  side) 
and  one  thermistor  on  the  silicon  anvil  to  which 
each  of  the  four  silicon  crystals  is  attached.  Vac- 
uum grease  is  used  to  insure  good  thermal  contact 
between  the  silicon  and  the  thermistors.  The  six 
thermistors  and  an  ultra-stable  calibrated  standard 
resistor  (2.5  kfi)  [22]  are  connected  in  series  and 
powered  with  a constant  current  source  ( = 1x10“^ 
A).  A precision  digital  voltmeter  is  used  to  read  the 
voltage  drop  across  each  thermistor  and  the  stan- 
dard resistor.  The  temperatures  are  read  every  4 s, 
so  the  temperature  for  a measurement  time  of  20  s 
is  an  average  of  five  measurements. 

Before  and  after  a critical  measurement,  the  six 
thermistors  are  placed  in  a constant  temperature 
bath  in  order  to  measure  offsets  between  the  ther- 
mistors. The  offsets  are  typically  a few  mK  and  are 
stable  to  less  than  1 mK  over  times  long  compared 
to  the  few  weeks  needed  for  a lattice  parameter 
measurement. 

The  temperature  of  the  laboratory  is  constant  to 
within  —0.1  K over  several  days.  For  the  measure- 
ment reported  here  the  spectrometer  was  isolated 
from  direct  air  currents  within  the  laboratory  by  a 
cloth  curtain.  The  curtain  has  subsequently  been 
replaced  by  an  insulated  chamber.  The  crystals  are 
covered  with  an  aluminum  thermal  shield  and  two 
20  mm  thick  styrofoam  thermal  shields.  The  x-ray 
tubes,  motors,  and  detectors  are  outside  the 
shields.  The  temperature  difference  between  the 
first  and  second  crystals  is  typically  0.1  K and  can 
be  varied  by  changing  the  temperature  of  the  water 
used  to  cool  the  x-ray  tubes.  The  temperature  dif- 
ferences between  the  second  crystals  is  only  a few 
mK.  Variation  of  the  temperature  differences  (the 
critical  quantity  for  relative  measurements)  over  a 
24  h period  is  typically  less  than  20  mK. 


separation  of  profiles  recorded  along  the  path  ( - ) 
and  path  ( -f- ).  The  profiles  are  measured  by  step- 
ping the  axis  through  N discrete  angles  (typically 
V = 110)  which  are  recorded  as  angle  interferome- 
ter fringe  numbers.  For  each  angle  the  number  of 
x-ray  photons  in  each  detector  is  counted  for  a 
fbced  time  (typically  10  s to  20  s).  The  time  per  scan 
is  40  min  to  80  min.  The  fringe  numbers  for  each 
profile  are  corrected  in  a point  by  point  manner  for 
the  temperature  difference  between  the  first  and 
second  crystals.  From  the  Bragg  equation  the  de- 
pendence of  the  diffraction  angle,  6,  on  tempera- 
ture, T,  is  easily  seen  to  be 

dd/dT  = —ai)tar\6  , 

where  ao  = linear  coefficient  of  thermal  expansion 
of  the  crystal.  For  silicon  00  = 2.56  x 10“ VK.  Let  9- 
and  9+  be  the  angular  settings  of  the  first  crystal 
for  the  path  ( — ) and  path  (-I-),  respectively,  the 
rotation  angle  reference  be  the  y axis,  and  clock- 
wise and  counter  clockwise  rotations  be  positive 
and  negative  respectively.  When  the  Bragg  condi- 
tion is  simultaneously  satisfied  at  crystals  1 and  2 
and  both  crystals  are  at  a reference  temperature. 
To,  (see  Fig.  la  and  Fig.  8), 

0-(r„)  = M7’o)-0B.(ro) 

and 

0+(7’o)  = — 0b2(T'o)  -I-  0bi(7o). 

When  crystal  1 is  at  temperature,  T,  and  crystal  2 
is  at  temperature,  T2,  then  these  two  equations  be- 
come 

An 

9-(To)  = [9s2(T2)  - 0Bi(ri)]  T2)  = 

An 

9-(TuT2)-^iT,-T,) 

and 

H R 

9,  (To)=-9^2iT2)  + 9Bi(T,)  + ^iT2-Ti)  = 

H f) 

9,{TuT2)  + ^iT2-T0, 


6.  Data  Analysis 

The  angular  offsets,  /3s  and  fin,  which  are  needed 
to  calculate  Ldjd  from  Eq.  (2)  are  the  angular 
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same  material. 
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The  angles  6-  for  the  path  ( — ) are  corrected  by 
— dO/dT  {T2  — T1)  while  the  angles  0+  for  the  path 
( + ) are  corrected  by  +dd/dT  (72-71).  Since  the 
profiles  are  recorded  in  fringe  numbers,  these  an- 
gular corrections  need  to  be  converted  to  fringes 
before  adding.  The  angle  measuring  interferometer 
is  positioned  within  a few  minutes  of  the  zero  de- 
gree angle  for  all  the  measured  offsets.  Using  the 
interferometer  equation  discussed  above,  the  angle 
to  fringe  conversion  factor  of  0.0401  sec/fringe  at 
zero  degree  angle  is  determined  and  used  to  con- 
vert the  angle  corrections  due  to  temperature  dif- 
ferences into  fringes. 

The  angular  offsets  between  the  profiles  on  the 
temperature  corrected  fringe  scale  have  been  de- 
termined by  two  different  methods.  In  the  first 
method  the  profiles  were  fit  with  a Lorentzian 
function  using  a nonlinear  least  squares  procedure 
in  which  the  position,  intensity,  width,  and  back- 
ground are  adjusted.  The  data  with  its  pronounced 
wiggles  and  central  spike  is  not  well  represented  by 
the  smooth  Lorentzian  function.  However,  by  tak- 
ing the  differences,  R±,  between  the  recorded  pro- 
files, I±,  and  the  smooth  fitting  functions,  L±,  the 
wiggles  of  the  two  profiles  are  emphasized. 

The  correlation  function,  C(A0),  for  the  two  sets 
of  residuals  is  computed. 

N 

C{A0)=^  R^id^)+R-(e,  + Ae), 

1 = 1 

where  the  sum  is  over  the  N data  points  of  the 
profile. 

The  value  of  A0  for  which  C(A0)  is  a maximum 
is  the  angular  separation  (offset  angle  (3)  of  the 
two  profiles.  The  uncertainty  in  the  measurement 
of  the  angular  separation  of  the  profiles  is  typically 
~6x  10“’  rad. 

In  the  second  method  the  dynamical  diffraction 
function  was  fit  to  the  recorded  profiles.  The  crys- 
tal thicknesses  were  fixed  at  0.455  mm  and  the 
structure  factors  were  taken  from  Refs.  [14  and 
15].  The  only  adjustable  parameters  in  the  fitting 
procedure  were  the  position,  the  intensity  and  the 
background.  The  angular  separation  between  the 
profiles  is  obtained  as  the  difference  between  the 
fitted  position  of  the  two  profiles.  The  measured 
angular  separations  obtained  with  the  two  methods 
agree  within  the  measurement  uncertainty.  Be- 
cause the  first  method  is  computationally  simpler 


and  more  easily  adapted  to  a small  computer,  it 
was  used  to  obtain  all  of  the  results  presented 
below. 

A single  lattice  comparison  run  usually  consists 
of  16  to  20  data  scans,  preceded  and  followed  by 
alignment  scans.  The  profiles  are  scanned  by 
repeating  the  following  sequence:  unknown  crys- 
tal—cw  rotation,  unknown  crystal  — ccw  rotation, 
standard  crystal— cw  rotation,  and  standard  crys- 
tal—ccw  rotation.  After  determining  the  angular 
separation  for  each  scan,  the  angular  separations  vs 
time  of  day  for  each  crystal  are  fit  with  a variable 
order  polynomial  (usually  3).  The  constant  (in 
time)  angular  offset  between  these  two  curves  is 
the  angular  offset  between  the  unknown  and  stan- 
dard crystals  for  one  data  set.  The  scatter  in  the 
offsets  is  typically  <6 x 10“’  rad  provided  the  crys- 
tals are  well  aligned  and  the  same  areas  of  the  crys- 
tals are  used  to  diffract  the  x-ray  beam. 

7.  Comparison  of  Two  Silicon  Samples 

Two  Si  samples  were  prepared  as  described  in 
Sec.  4 from  material  acquired  by  PTB  and  NIST 
for  absolute  lattice  parameter  measurements.  The 
NIST  sample  was  supplied  by  Dow  Chemical  and 
was  a slab  adjacent  to  the  NIST  x-ray/optical  inter- 
ferometer. The  PTB  material  was  supplied  by 
Wacker-Chemitronic  and  labeled  WASO  17  by 
PTB.  The  relative  difference  between  WASO  17 
and  WASO  4.2  (the  PTB  x-ray/optical  interferome- 
ter crystal)  has  been  measured  to  be  2.5±1  x 10““ 
[23](WASO  17 -WASO  4.2).  The  first  long  crystal 
was  prepared  from  material  supplied  by  Monsanto. 

These  two  samples  were  chosen  because  the  lat- 
tice spacing  of  the  NIST  sample  has  been  mea- 
sured absolutely  and  the  lattice  spacing  of  WASO 
17  was  measured  relative  to  the  WASO  4.2  sample 
which  was  measured  absolutely  at  the  PTB.  In  the 
case  of  the  NIST  sample,  further  refinement  of  the 
absolute  lattice  parameter  measurement  is  continu- 
ing [24,  25].  Because  the  absolute  lattice  measure- 
ments have  not  yet  achieved  0.01  ppm  uncertainty, 
they  do  not  provide  as  definitive  a test  of  the  Ad 
measurements  as  is  desirable.  The  published  re- 
sults of  the  absolute  lattice  spacing  measurements 
are 


^ Certain  commercial  equipment,  instruments,  or  materials  are 
identified  in  this  paper  to  specify  adequately  the  experimental 
procedure.  Such  identification  does  not  imply  recommendation 
or  endorsement  by  the  National  Institute  of  Standards  and 
Technology,  nor  does  it  imply  that  the  materials  or  equipment 
identified  are  necessarily  the  best  available  for  the  purpose. 
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^^22o(WASO  4.2)  at  22.5  °C  = 0.192015560  (12)  nm 
(±6.2x10-*) 

^/22o(NIST)  at  22.5  ‘’C  = 0.19201554  (2)  nm 
(±10.4x10-*) 

which  lead  to 

Ac//(/(NIST-WASO  4.2)  = - Ixl0-’±1.2xl0-^ 

In  order  to  develop  confidence  in  the  Aif  mea- 
surements, we  have  measured  the  long  term  repro- 
ducibility by  recording  an  extensive  set  of  data  over 
more  than  one  half  year  during  which  the  crystals 
were  realigned  on  the  spectrometer  several  times. 
In  addition,  measurements  were  taken  using  two 
different  wavelengths  (Ag  Ka  and  Mo  Ka  x-ray 
radiation)  and  with  the  long  first  crystal  in  the  two 
possible  orientations  (the  0°  and  the  180°  orienta- 
tion). 

In  Fig.  7 the  measured  values  of  ^dld  are  plot- 
ted vs  time  with  different  symbols  for  the  different 
sources  and  different  first  crystal  orientations. 
Note  the  long  time  span  over  which  the  Ag  Ka,  0° 
orientation  measurements  were  taken.  In  Table  1 
the  average  numerical  values  for  the  different 
sources  and  first  crystal  orientations  are  presented. 
From  both  the  table  and  the  figure,  the  magnitude 
of  systematic  effects  related  to  the  source  and  first 
crystal  orientation  is  estimated  to  be  —7x10-'^. 
Table  2 provides  a summary  of  uncertainty  contri- 
butions resulting  from  systematic  effects  associated 
with  1)  the  wavelength  and  first  crystal  orientation, 
2)  crystal  temperature  measurements,  3)  crystal 
misalignments,  and  4)  periodic  nonlinearity  in  the 
Michelson  angle  interferometer. 
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Fig.  7.  Comparison  of  NIST  and  WASO  17  crystals;  + AgKa, 
first  crystal  0°;  □ AgKa,  first  crystal  180°;  O MoKa,  first  crystal 
0°;  A MoKa,  first  crystal  180°. 


Table  1.  Comparison  of  NIST  and  WASO  17  samples 


Source 

First  crystal 
orientation 

Adid  X lO" 

NIST-WASO  17 

AgKa 

0° 

0.508  ±0.379 

Ag  Ka 

180° 

1.236  ±0.393 

Mo  Ka 

0° 

1.050  ±0.262 

Mo  Ka 

180° 

1.251  ±0.310 

Table  2.  Estimated  uncertainties  from  systematic  effects 

Systematic  effect 

Contribution  to  AdId  x 10** 

A and  1st  crystal  orientation 

0.7 

Temperature 

0.3 

Crystal  misalignment 

0.2 

Interferometry  nonlinearity 

0.5 

The  uncertainty  in  the  measurement  of  the  crys- 
tal temperature  differences  is  less  than  1 mK.  This 
value  was  determined  by  periodically  placing  all  of 
the  thermistors  in  a constant  temperature  bath  as 
explained  in  Sec.  5.  A 1 mK  uncertainty  in  the  crys- 
tal temperature  differences  leads  to  an  uncertainty 
of  3 X 10"’  in  ^dld.  The  uncertainty  introduced  by 
crystal  misalignment  is  discussed  in  detail  in  Ap- 
pends A and  is  estimated  to  be  less  than  2x  10"'^. 

Heterodyne  Michelson  interferometers  are  prone 
to  sub-periodic  nonlinearity  resulting  from  imper- 
fect separation  of  the  two  frequencies  by  the  pri- 
mary beam  splitter  [26, 27].  The  periodic 
nonlinearity  can  be  evaluated  by  pressure  scanning 
the  interferometer  and  recording  the  fringe  ad- 
vance vs  pressure  increase.  The  maximum  ampli- 
tude of  the  nonlinearity  is  estimated  to  be  4 x 10“’ 
rad.  This  means  that  the  angular  separation  of  two 
points  on  a diffraction  profile  which  are  separated 
by  0.5  interferometer  fringe  (4  x 10"’  rad)  might  be 
in  error  by  8 x lO"''  rad.  However,  because  the  pro- 
files are  typically  25  interferometer  fringes  wide, 
the  influence  of  the  interferometer  periodic  error 
on  the  profile  peak  position  is  significantly  reduced 
by  averaging.  In  addition,  the  phase  relation  of  the 
x-ray  profile  to  the  angle  interferometer  fringes 
changes  so  that  uncertainty  due  to  periodic  nonlin- 
earity becomes  part  of  the  statistical  uncertainty. 
We  estimate  that  the  uncertainty  contribution  to 
Ldld  due  to  the  systematic  effect  associated  with 
the  periodic  nonlinearity  of  the  angle  interferome- 
ters is  not  larger  than  5 x 10"''.  The  angle  interfer- 
ometer is  being  modified  in  order  to  reduce  the 
periodic  nonlinearity  by  approximately  a factor  of  5. 
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A final  value  for  the  NIST-WASO  17  compari 
son  was  obtained  by  considering  each  of  the  entries 
in  Table  1 to  be  an  independent  measurement  of 
equal  weight.  The  statistical  uncertainty  was  com- 
bined with  the  uncertainties  from  systematic  effects 
in  Table  2 to  obtain  a value  for  ^dld  (NIST- 
WASO  17)  = (1.037  ± 1.0)  X 10-“*.  By  combining  this 
value  with  the  WASO  17-WASO  4.2  difference 
noted  above,  one  obtains 

^d|d  (NIST-WASO  4.2)  = (3.5  ± 1.4)  x lO'*. 

This  value  is  only  slightly  outside  the  1 a uncer- 
tainty of  the  absolute  measurements.  In  addition,  it 
should  be  remembered  that  the  NIST  absolute 
value  is  a preliminary  result. 

8.  Conclusions 

A lattice  comparison  facility  has  been  estab- 
lished at  NIST  which  is  capable  of  measuring  crys- 
tal lattice  spacing  difference  with  an  uncertainty 
<0.01  ppm.  By  using  crystals  of  equal  thickness, 
the  recorded  profiles  exhibit  fine  structure  which 
permits  more  precise  measurement  of  the  small  an- 
gular offsets  between  profiles.  The  spectrometer  is 
designed  to  permit  easy  interchange  of  crystal  sam- 
ples and  the  comparison  of  four  samples  in  one 
setup.  Procedures  are  provided  for  precise  align- 
ment of  the  crystals  so  that  alignment  errors  con- 
tribute <2x  10“^  A comparison  of  samples  whose 
lattice  spacings  have  been  absolutely  measured  is 


consistent  with  the  absolute  measurements.  In  the 
near  future  other  crystal  samples  (including  Ge) 
destined  for  x-  and  gamma-ray  diffraction  will  be 
compared  using  the  spectrometer  and  the  tech- 
nique described  here. 

9.  Appendix  A.  Derivation  of  the 
Equation  for  Calculation  of  A</ 

If  the  lattice  planes  in  crystal  1 and  2 are  all 
perfectly  parallel  with  the  axis  of  rotation  and  the 
x-ray  beams  lie  in  a plane  normal  to  the  axis  of 
rotation,  then  the  difference  in  lattice  spacing  be- 
tween crystal  1 and  2 can  be  derived  from  Bragg’s 
law  (see  Sec.  2).  In  reality,  however,  the  alignment 
of  the  crystals  and/or  the  x-ray  beams  is  never  per- 
fect. In  this  appendix,  equations  for  the  calculation 
of  AJ  which  include  misalignment  terms  are 
derived.  The  x-ray  beams  are  described  by 
wavevectors,  k,  and  the  crystals  are  described  by 
reciprocal  vectors,  t.  The  derivation  is  based  on 
the  fact  that  the  angular  positions,  0+  and  0-,  of 
the  symmetric  rocking  curve  profiles,  measured 
along  path  (-I-)  and  path  ( — ) in  the  two  crystal 
spectrometer,  are  equal  to  the  two  particular  rota- 
tion angles  0,  for  which  the  Laue-conditions  for 
path  (-f)  and  path  (-),  are  simultaneously  satis- 
fied at  crystal  1 and  2 (see  Figs.  8 and  9). 

Larson  [28]  has  made  a similar  geometrical  anal- 
ysis of  a lattice  parameter  comparison  using  Bragg 
geometry.  We  first  derive  the  equation  for  l^d/d  for 
the  ideal  case  involving  perfect  crystals  and  no  mis- 
alignment of  the  crystals  or  the  x-ray  beams.  The 


Fig.  8.  Geometr)'  of  the  lattice  parameter  measurement  with  no  crystal  mis- 
alignment. 
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z 


Fig.  9.  Geometry  of  the  lattice  parameter  measurement  including  crystal 
misalignment. 


more  experienced  reader  may  choose  to  skip  over 
this  section.  Next,  we  derive  the  equations  for  the 
nonideal  case  and  apply  them  to  “alignment  scans” 
and  “data  scans.” 

9.1  Ideal  Case 

First  consider  the  ideal  case  of  no  misalignment 
of  the  x-ray  beams  or  tilting  of  the  crystals.  The 
geometry  of  the  measurement  is  then  entirely  con- 
fined to  the  xy  plane  and  is  as  shown  in  Fig.  8.  The 
coordinate  system  is  such  that  the  z-axis  is  parallel 
to  the  0-axis  (i.e.,  the  axis  of  the  rotating  table), 
the  X -axis  is  parallel  to  the  lattice  planes  of  interest 
in  crystal  2,  the  y-axis  is  normal  to  the  lattice 
planes  in  crystal  2,  and  the  normal  to  the  lattice 
planes  in  crystal  1 lies  in  thexy  plane.  The  rotation 
angle  reference  is  taken  to  be  the  y-axis  and  the 
0’s,  the  angles  of  rotation  of  the  first  crystal,  are 
positive  for  clockwise  rotation  and  negative  for 
counterclockwise  rotation. 

The  x-ray  beams  are  described  by  their  wave  vec- 
tors: k-,  for  the  incident  rays;  kl,  for  the 
singly  diffracted  rays;  and  k-,  kl  for  the  double 
diffracted  rays.  Using  the  angles  shown  in  Fig.  8, 
the  doubly  diffracted  rays  can  be  written  as 

It"  =27r/A(cosa-,  sina",  0), 

k- =2Trl\{cosal,  sina+,  0). 

The  crystals  are  described  by  their  reciprocal  lat- 
tice vectors  ti(0)  for  crystal  1 and  t2  for  crystal  2 
which  can  be  written  as 


Ti(0)=^  (sin0,  COS0,  0), 

= ^ (0,1,0). 

The  planes  at  the  two  ends  of  crystal  1 are  assumed 
to  be  parallel  and  to  have  the  same  lattice  spacing. 
When  the  Laue  condition  is  satisfied  at  crystal  1 

kL  —k-  = - Ti  (0-)  and  Ar+  -*+  = T((0+).  (3) 

Similarly  when  the  Laue  condition  is  satisfied  at 
crystal  2 

k--k'-  = ndindkl-k'^=-Ti.  (4) 

Since  the  wave  vectors  have  equal  magnitude, 
2-77/4,  Eqs.  (3  and  4)  imply 

an<lT,(e.)-(i5+t;)  = 

||T,{e.)|’.  (5) 

By  substituting  the  reciprocal  lattice  vectors  and 
the  wave  vector  with  a"  = a"  = 0b2  = Bragg  angle  at 
the  crystal  2 into  Eq.  (5)  and  using  the  Bragg  equa- 
tion, A = 2d2  sin02,  to  eliminate  the  wavelength,  two 
equations  for  the  lattice  spacing  difference  are 
obtained. 
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— sin0-cot0B^  - (1 -cos0-)~ 
d\ 

— 0-COt0B2-^  0- 

= -sin0+cot0B2-(l  -cos6+)~ 
di 

0+COt0B2-|  9l 

where  terms  up  to  second  order  in  0^  have  been 
retained.  These  two  equations  can  be  added  to  ob- 
tain 


d2-di 

di 


cot0B2-^(0--0+)^  (6) 


where  we  have  used  the  fact  that  d+  — — 6-. 

For  comparisons  in  which  the  first  and  second 
crystal  are  made  out  of  the  same  material 
[(^/2 < 10“^],  the  last  term  is  negligible  as 
the  following  inequality  shows 


d2-dx 

2 1 

d2-dy 

d, 

^2 

dy 

tan^0B2== 


I if  0b2<54°. 

If  crystal  2 is  alternately  the  standard  and  unknown 
crystal,  the  Eq.  (6)  becomes 

~ cot0Bs  where  /3s  = 0-s - 0+s 

~ cot0Bu  where  /3u=  0-u-  0+u 

Ul  / 

These  two  equations  immediately  lead  to  Eq.  (2). 


with  respect  to  the  plane  normal  to  the  rotation 
axis  introduce  correction  terms  to  the  measured 
lattice  spacing  difference.  In  the  analysis  which 
takes  into  account  misalignments  and  crystal  im- 
perfections, the  geometry  of  Fig.  9 is  used. 

The  coordinate  system  is  identical  to  that  used 
for  the  ideal  case,  i.e.,  the  z axis  is  parallel  to  the 
axis  of  rotation  and  thex  axis  parallel  to  the  lattice 
planes  in  crystal  2.  However,  the  wave  and  recipro- 
cal lattice  vectors  no  longer  lie  in  the  xy  plane.  The 
wave  vector  of  the  double  diffracted  x-rays  can  be 
put  in  the  following  form: 


k"  = — 

- A 


(cosa±  cos5±,  -t-sinaj;,  cosa±  sin5±),  (8) 


where  8±  is  the  angle  between  thej^y  plane  and  the 
projection  of  the  wavevector  it±  on  thexz  plane. 


The  reciprocal  lattice  vectors  associated  with  path 
(-F)  and  path  (-)  in  crystal  1,  t+  and  t-,  are  no 
longer  parallel  and  of  equal  magnitude.  The  two 
vectors  t±  can  be  described  by  the  six  parameters 
d±,  p±  and  6±  as  follows: 

T±  (sin(0  - (/>3:)  cosp±, 

cos(0  - <^±)  cosp±,  sinp±),  (9) 

where  d±  is  the  lattice  spacing  in  the  diffracting  re- 
gion of  crystal  1 corresponding  to  path  ( ± ),  p±  is 
the  angle  between  t±  and  thexy-plane,  and  is 
the  particular  angle  of  rotation  6 for  which  t*  is 
parallel  to  the  yz  plane.  For  a perfect  first  crystal, 
d+  =d-,  p+  =p-,  and  </>+  =4)-.  In  crystal  2 the  re- 
ciprocal lattice  vector,  t5,  is  the  same  for  both  paths 
and  is  given  by 

t2  = ^ (0,  cosp,  sinp),  (10) 

«2 


^s^^_A^COt0Bs  (7) 

Ul  2 


9.2  Non-Ideal  Case 

Misalignment  of  the  planes  of  crystals  1 and  2 
with  respect  to  the  rotation  axis,  imperfections  in 
the  first  crystal  (i.e.,  lattice  spacing  gradient  and 
nonparallelism  of  the  planes  at  the  two  ends  of  the 
crystal),  and  the  misalignment  of  the  x-ray  beams 


where  di  is  the  lattice  spacing  in  the  diffracting  re- 
gion of  crystal  2,  and  p is  the  angle  between  n and 
xy-plane.  Note  that  tz  is  in  theyz  plane. 

For  path  (±),  the  relations  between  the  wave 
vectors  k±  of  the  incident  rays,  kl  of  the  single  dif- 
fracted rays,  and  k±  of  the  double  diffracted  rays 
are  given  by  the  Laue  conditions: 


k±-k±=±T±  (crystal  1) 

(11) 

k'i—kl= +72  (crystal  2) 

(12) 
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Since  the  wave  vectors  have  equal  magnitude,  Now  by  inserting  Eqs.  (16),  (17),  (18),  and  (19)  into 
Ztt/A,  Eqs.  (11)  and  (12)  imply  that  Eq.  (14)  the  following  equation  can  be  derived: 


Tz  -/cS  = +1  Irzp 


r^'{Ti±kl)=^T±f. 


(13) 


di-d^ 


= ±(0±  -<f>±)  cot0B+e±  , 


(20) 


(14) 


where  e±  is  the  second  order  correction  term 


When  the  Laue  condition  is  satisfied  at  crystal  2, 
the  angles,  a±,  can  be  expressed  as  a±  = 0B2  + e± 
where  e±  are  the  small  angular  corrections  to  the 
second  crystal  Bragg  angle  because  of  crystal  imper- 
fections and  misalignments.  Thus  Eq.  (8)  becomes 

k±=^  (cos(0b2+ e±)  cos5±, 

A 

+ sin(0B2  + e±),  cos(0B2  + e±)  sin5±).  (15) 

In  order  to  establish  a relation  between  the  lat- 
tice spacing  d of  crystal  2 and  the  rotation  angles 
9 = 0±,  for  which  the  Laue-conditions  correspond- 
ing to  path  ( ± ) are  simultaneously  satisfied  at  crys- 
tal 1 and  crystal  2,  Eqs.  (9),  (10),  and  (15)  are 
expanded  to  second  order  in  the  angles  6±  — (f)±,  p±, 
p and  and  to  first  order  in  the  correction  angle 
€±: 

t.=22  (0.-.^.,  pl,p.)(16) 


e±=  - j (9. (p.  -p)  (p.-3p) 

±5±  (pi -p)  cot0B2  . (21) 

The  angle  8±  in  Eq.  (21)  describes  the  inclination  of 
the  double  diffracted  beam  relative  to  the  Ay -plane, 
(see  Fig.  9).  If  crystal  1 and  crystal  2 are  perfectly 
parallel,  the  wavevectors  k±,k'±,  and  k'i  are  copla- 
nar  and  the  inclinations  5±  are  determined  simply 
by  the  positions  of  the  source  and  detector  slits. 
When  the  two  crystals  are  not  perfectly  parallel,  8± 
are  also  dependent  on  the  misalignment  angles  p± 
and  p. 

Let  £,  t,  and  €"  be  projections  along  thex-axis  of 
thex-ray  paths  from  the  source  slit  to  crystal  1,  from 
crystal  1 to  crystal  2,  and  from  crystal  2 to  the  detec- 
tor slit,  respectively.  The  distances  are  shown  in  Fig. 
1 and  for  our  spectrometer  are: 

€ =(25.4  COS0B  + 8.25)  cm 

€'  =8.25  cm 

f'  = (28.5  cos0b)  cm 


(17) 


((1 5|)  cot0B2-ei:,  +l  + e± 


COt0B2, 


If  AZ±  denotes  the  difference  in  Z- coordinate  be- 
tween the  detector  slit  and  source  slit  of  path  ( ± ), 
then  the  inclination  A*  relative  to  the  jy -plane  of 
the  straight  line  joining  the  two  slits  is  given  by: 


A^  = 


AZ^ 

€ + £'  + €"  • 


(22) 


6±cot0B2)-  (18) 


In  the  last  equation,  the  wavelength  A has  been 
eliminated  using  A =2^/2  sin 0b2.  By  insertion  of  Eqs. 
(17)  and  (18)  into  Eq.  (13),  an  expression  for  the 
correction  angles  e±  is  easily  obtained: 


The  change  in  Z -coordinate  as  one  moves  from  the 
source  slit  to  the  detector  slit  along  path  ( ± ) can  be 
expressed  by  combining  the  Z -components  of  the 
wave  vectors  k±,  k±,  and  k'i,  normalized  to  unity 
and  the  distances  traveled: 

AZi  =^  [(it.)z^+ (it:).r  + (it:)z€" 


€±=2  p^  tan0B±p5i 


(19) 


1 

COS0B  " 


(23) 
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From  Eqs.  (11),  (12),  and  (15)  and  using  Eqs.  (9) 
and  (10),  the  Z -components  of  the  wave  vectors 
are  found  to  second  order  in  p±,  p,  and  5±: 


^(it±)z  = 5±  cos0b±2  (p-p±)  sin0B  (24) 
^ (ifci)?  = 5±  cos0B±2p  sin0B  (25) 

Ztt 

^(it:),  = 5^  COS0B.  (26) 


By  insertion  of  Eqs.  (22),  (24),  (25),  and  (26)  into 
Eq.  (23)  a relation  between  and  A±  is  obtained: 

8^=A^±2  + + P^- ( + P)  • 

Using  this  expression,  the  unknown  angles  S±  can 
be  eliminated  from  Eq.  (21): 

e±=  (0--(^i)^  + 2 ip±-p)  {aip±+bip)± 


A±(p±- p)  cotOs  (27) 


where 

e +£'+("  4 


e+e'  + t"  4 


Even  though  the  distances  € and  ("  depend  on  the 
Bragg  angle  0b,  the  numerical  values  for  ai  and  bi 
are  essentially  identical  for  the  two  x-ray  sources 
(0b  = 16.9°  for  Ag  Ktti  and  0b  = 21.6°  for  Mo  Koi). 

When  crystal  1 and  2 are  both  of  the  same  mate- 
rial, an  upper  limit  of  the  first  term  in  Eq.  (27)  is 
easily  estimated: 


5 


d-d^  2 
d^ 


<io-‘° 


io-“>> 


d-d^ 

2 1 

1 

1+ 

d^ 

""2 

tan^0B~|(0.  - <t>±y, 


if  0b<54°. 


This  term  is  two  orders  of  magnitude  smaller  than 
the  intended  uncertainty  (a  few  parts  in  10“)  in  the 
lattice  spacing  comparison  and  is  therefore  ne- 
glected in  the  following.  The  remaining  two  terms 
in  Eq.  (27)  are  both  zero,  if  the  lattice  planes  in 
crystal  1 and  2 are  parallel,  i.e.,  p+-p  = 0 and 
p_-p=0.  The  aim  of  the  alignment  scan  is  to 
measure  p+-p  and  p_—p  and  then  to  adjust  the 
tip  of  crystal  2 so  as  to  make  p^—  p and  p_  -p  as 
close  to  zero  as  possible. 

9.2.1  The  Alignment  Scan  In  the  alignment 
scan,  only  one  path,  path  ( + ) or  path  ( - ),  is  used. 
Two  rocking  curves  are  measured  simultaneously: 
one  curve  using  the  upper  (up)  position  and  one 
using  the  lower  (lo)  position  of  the  detector  slit. 
For  the  path  (-),  the  observed  difference  A0-  in 
the  centers  of  the  two  rocking  curves 
(A0- = 0-up- 0-io),  the  change  in  the  inclination 
(A-up-A_io)  when  the  detector  slit  is  moved  from 
the  upper  to  the  lower  position,  and  the  difference 
in  misalignment  of  the  two  crystals,  p-  -p,  are  re- 
lated. The  following  equation  is  derived  directly 
from  Eqs.  (20)  and  (27): 


(28) 


Similarly,  for  path  ( + ) it  follows  that 


P+- 


(29) 


(A-up-A-io)  and  (A+up-A+io)  can  be  estimated 
from  mechanical  measurements  of  the  spectrome- 
ter and  are  ~9  x 10"^  rad.  Due  to  bending  of  crys- 
tal 1,  p+  and  p-  are  not  equal.  Consequently,  the 
sample  crystal  can  only  be  aligned  for  one  of  the 
two  paths.  If  the  tilt  of  crystal  2 is  adjusted  so  that 
|A0-|  =2x10”“  rad,  then  typical  values  for  |A0+| 
are  1 x 10“’  rad.  These  values  imply  that 

|p--p|  = 2xl0““ 

|p+  — p|  = l X 10'^  . 

The  magnitude  of  the  individual  p’s  can  also  be 
controlled  and  estimated.  By  rotating  the  1st  crys- 
tal by  180°,  the  tilt  of  the  first  crystal  can  be  ad- 
justed so  that  p+  and  p-  are  < 1 x 10“^  rad. 

9.2.2  The  Data-Taking  Scan  In  the  data- 
taking  scan,  two  rocking  curves  are  simultaneously 
measured  with  the  detector  slits  fixed  in  the  middle 
position:  one  curve  using  path  (-F),  and  one  using 
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path  ( - ).  We  define  the  average  lattice  spacing,  di, 
in  crystal  1 as 

, 2d^d- 

and  average  Eqs.  (20)  and  (27)  over  path  ( + ) and 
path  ( - ) to  obtain 

(0.-0-  + Tjo)  cot 0B  + e,  (30) 

where 

^ ={(p+  -p)(«ip+  +^ip)  + (p-  -p)(«iP-  +^ip)}  + 

I {^+(p+ -p)-^-(p- -p)}  cot0B  (31) 

r7o=</._-<^..  (32) 

The  angle  rjo  describes  a bending  of  crystal  1 around 
the  0-axis  and  is  independent  of  the  particular  sam- 
ple used  as  crystal  2.  Obviously,  the  appearance  of 
the  unknown  angle  770  in  Eq.  (30)  excludes  the  pos- 
sibility of  using  crystal  1 as  an  absolute  reference 
crystal. 

The  two  crystals  being  compared  are  placed  on 
the  translation  table  and  are  used  alternately  as 
crystal  2.  One  of  the  crystals  serves  as  the  standard 
(s)  and  one  as  the  unknown  (u).  The  difference  in 
lattice  spacing,  ^d  =d^  — dn  for  the  two  sample  crys- 
tals can  be  derived  from  Eqs.  (30)  and  (31). 

cot0B  + Ae  (33) 

where: 

/3s=  0-s  - 0+s 

/3u=  0-U-0  + U 

Ac  =[^2(0+  +P-)  + ^2(Ps  + Pu)  — 2 

(A+ - A-)cot0B]  (ps-pu)  = 

[a2(p+  +p-)  + 2Z22Pu-|  (A+  - A_)cot0B] 

(Ps  - Pu)  + bi  (ps  - pu)^  (34) 

where  a2  = bi-ai=  -0.08  and  b2=  -2bi=  -0.30. 

In  order  to  make  Ae  small,  the  tilts  of  the  crystals 
are  adjusted  so  that  p’s  are  <1x10“’  (see  Sec. 


9.2.1).  Care  is  taken  to  center  the  source  and  detec- 
tor slits  on  the  plane  of  dispersion  in  order  to  make 
A+-A-  small,  but  mechanical  measurements  sug- 
gest that  the  magnitude  of  the  A’s  is  ^ 3 x 10"“*.  The 
correction  term  can  be  estimated  by  measuring 
()8s  - /Su)  as  a function  of  (p,  - pu)  where  pu  is  kept 
fixed  and  ps  is  varied.  For  ps-pu  ^5x10“’,  Eqs. 
(33)  and  (34)  can  be  approximated  by 

^=-(^^)cotft,  + C(p.-p.) 

where 

C =fl2(p+  + P-)  + 2b2Pu~2  (^+  — A-)cot0B. 

A plot  of  (73s  - /3u)  vs  (ps  - Pu)  is  shown  in  Fig.  10  and 
has  an  approximate  slope  of  -3.6x  10“'*  in  reason- 
able agreement  with  the  estimates  of  the  A’s  from 
the  mechanical  measurements  and  the  estimates  of 
the  individual  p’s. 


(ps-pu)  X 10®  * radians 


Fig.  10.  Measured  data  profile  offsets,  /?,-  /3u,  as  a function  of 
the  measured  crystal  tilts  (ps-pu).  See  text  for  more  explanation. 

Before  and  after  data  scans  are  recorded  using  a 
standard  and  unknown  crystal,  alignment  scans  are 
recorded  for  the  standard  and  unknown  crystal  also. 
The  quantity  ps-pu  is  determined  from  the  align- 
ment scans  and  is  maintained  at  a value  <4x  10“* 
rad.  Thus  the  correction  term  Ae  < 2 x 10“’.  Correc- 
tion of  the  data  for  measured  values  of  Ae  does  not 
reduce  the  statistical  spread  of  repeated  measure- 
ments. We  have  thus  chosen  to  increase  the  relative 
uncertainty  of  the  final  results  by  2x  10“’,  but  not 
to  correct  individual  measurements.  Although  the 
correction  term  is  small,  it  could  be  reduced  by  an 
order  of  magnitude  by  making  the  source  and 
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detector  slits  adjustable  about  the  plane  of  disper- 
sion. Modifications  to  the  spectrometer  to  permit 
this  adjustment  are  in  progress. 

10.  Appendix  B.  Electronics  Associated 
with  the  Interferometer 

The  layout  of  the  Michelson  angle  interferome- 
ters is  shown  in  Fig.  2 and  is  briefly  described  in 
Sec.  3.  In  this  appendix  the  electronics  used  to 
measure  the  whole  and  fractional  fringes  and  to 
servo-control  the  rotation  axes  are  described.  A 
block  diagram  of  the  electronics  is  provided  in  Fig. 

11. 

The  1.8  MHz  TTL  level  signals  from  the  laser 
and  the  detector  are  first  converted  to  low 
impedance  signals  by  buffer  amplifiers.  The  laser 
beat  signal  is  then  passed  through  a 2-way  90°  split- 
ter with  0°  and  90°  split  yielding  sin  {(ot)  and  cos 
{(lit),  respectively.  The  detector  beat  signal  is 
passed  through  a 2-way  0°  splitter  yielding  two  sin 
{iot  + dj))  signals.  The  laser’s  cos  (<ur)  and  one  of 
the  detector’s  sin  {o)t  + 4>)  are  combined  in  a mixer 
which  gives  signals  at  the  sum  and  difference  fre- 
quencies. The  difference  frequency  signal  is  sin  ({) 
and  the  sum  frequency  signal  is  sin  (2cot  -I-  </>)  which 
is  attenuated  by  a low  pass  filter.  Meanwhile  the 
sin  {(ot)  from  the  laser  is  combined  with  the  other 
sin  (cot  + (f))  of  the  detector  in  another  mixer  which 
after  filtering  gives  cos  cj).  The  sin  ^ and  cos  4)  sig- 
nals are  next  amplified  to  restore  their  signal 
strengths  which  have  deteriorated  due  to  losses  in 
the  splitters,  mixers,  and  filters. 

Ninety  degree  changes  in  4)  which  correspond  to 
— 0.04  second  rotation  of  the  rotary  table  are  de- 
tected by  sending  the  sin  4>  and  cos  4>  signals  to  an 
up/down  counter.  Smaller  changes  in  4>  are  mea- 
sured by  comparing  4)  to  a standard  angular  signal 


(sin  6 and  cos  6)  obtained  as  the  output  of  a digital 
vector  generator.  A vector  generator  with  a few 
0.1°  sensitivity  provides  a total  angular  sensitivity  of 
a few  X 10"''  seconds  with  the  vector  generator  ad- 
dress being  directly  related  to  the  fringe  fraction. 
By  taking  the  difference  of  the  product  of  sin  4>  and 
cos  6 and  the  product  of  cos  </>  and  sin  6 the  signal 
sin(4>  - 0)  is  obtained.  This  is  an  appropriate  error 
signal  for  locking  the  angle  of  the  rotary  table  via  a 
piezoelectric  transducer  to  the  position  where 

(f>  = e. 
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